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ABSTRACT 

Surface Acoustic Wave (SAW) technology is low cost, 
rugged, lightweight, extremely low power and can be used 
to develop passive wireless sensors. For these reasons, 
NASA is investigating the use of SAW technology for 
Integrated Vehicle Health Monitoring (IVHM) of aerospace 
structures. To facilitate rapid prototyping of passive SAW 
sensors for aerospace applications, SAW models have been 
developed. 

This paper reports on the comparison of three methods 
of modeling SAWs. The three models are the Impulse 
Response Method a first order model, and two second order 
matrix methods; the conventional matrix approach, and a 
modified matrix approach that is extended to include 
internal finger reflections. The second order models are 
based upon matrices that were originally developed for 
analyzing microwave circuits using transmission line 
theory. Results from the models are presented with 
measured data from devices. 

Keywords'. Surface Acoustic Wave, SAW, transmission line 
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1 INTRODUCTION 

First order models of SAW devices are based upon the 
Impulse Response [1, 2]. These models do not take into 
account second order effects such as internal reflections, 
frequency shifts, or allow for any physical arrangement 
other than equal electrode widths and spaces. For more 
accurate results, a matrix based approach was developed 
[3]. This approach has been further refined and modified to 
include internal finger reflections [4]. The reflections occur 
when the thickness of the metallization is sizeable enough 
to result in significant reflections. The extensions are based 
upon matrices that were originally developed for analyzing 
microwave circuits using transmission line theory. The 
modifications are accomplished by breaking up the SAW 
device into zones, where the area under a metalized region 
is treated as one zone, and the area without metallization is 
treated as another zone. The impedance discontinuities that 
occur at the edges of the metal fingers enable the simulation 
of the internal reflections of the mechanical acoustic wave. 
The modifications also enable incorporation of the different 
velocities for each region, which produces a more accurate 
characterization of the frequency response of the device. 


In this paper, the three modeling methods are briefly 
discussed, then a comparison of the three models with 
results from prototype devices are presented. 

2 FIRST ORDER MODELING USING THE 
IMPULSE RESPONSE METHOD 

The Impulse Response method [1] was used as the 
baseline for modeling the SAW device. This method is 
valid only for transducers where at least one of the two 
Inter-Digitated Transducers (IDTs) has a constant aperture 
or finger overlap [3]. This first order model includes both 
the mechanical and electrical behavior of SAW devices. It 
calculates the frequency response, the loss of the system, 
the admittance, and parameters for circuit simulators. This 
model assumes constant and equal spacing and finger 
widths. A simple circuit model (Fig. 1) can be used to 
convey the basic elements of the Impulse Response Model. 
The figure shows the source voltage and both the source 
and load impedances which are not part of the model. In 
the circuit model Ct is the total capacitance, B a (f) is the 
acoustic susceptance and G a (f) is the radiation conductance. 
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Fig. 1. (a) Basic SAW delay line and (b) the circuit model 
used in the Impulse Response Modeling. Ct is the total 
capacitance, B a (f) is the acoustic susceptance, and G a (f) is 
the radiation conductance. 


The frequency response of a SAW device can be 
calculated by using the Impulse Response Model and is 
given by 
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Where k is the piezoelectric coupling coefficient, C s is 
the capacitance per finger pair and unit length, H ;l is the 
aperture or overlap height of the fingers, fo is the center or 
synchronous frequency, N p is the number of finger pairs, f 
is the frequency, D is the delay length between the IDTs, 
and X is given by: 


v Ar M-fo) 

X = N p n . ( 2 ) 

Jo 

3 CONVENTIONAL MATRIX METHOD 

The matrices discussed in this work were originally 
developed for analyzing microwave circuits using 
transmission line theory. The methodology utilizing 
transmission matrices was modeled on the approach given 
by Campbell [3]. This method is based upon the Mason 
equivalent circuit using the crossed field technique (Fig. 2). 
Where for modeling purposes an IDT can be modeled as a 
single entity with an electrical port, and two acoustic ports. 
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Fig. 2. Transmission matrix model of an IDT. 

This allows the acoustic waves (W;) and electrical 
parameters (a; and b;) to be related through the use of 
transmission matrix T in: 
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The transmission matrix is in turn broken up into sub- 
elements, given by 
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Given the T matrix for an IDT, calculations for a SAW 
delay line or filter can be performed. The matrix for a 
SAW delay line is simply the multiplication of a 4x4 sub- 
matrix (elements tn, t [2 , t 2 i, and t 22 ) for the two IDTs and a 
matrix for the delay in between (Fig. 3). 
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Fig. 3. Transmission matrix model of a complete SAW 
delay line comprising of two IDTs and the delay between. 

The delay matrix is modeled after an acoustic 
transmission line as well. The delay matrix is given by: 
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where X is the wavelength at the synchronous frequency 
and d is the distance between the reference planes, or in this 
case the center of the two IDTs. Therefore complete SAW 
device matrix is given by: 

[SAW(f)]=[T l {f)D l {f)T 2 {f)~\. (6) 

4 MODIFIED MATRIX METHOD 

For more accurate results, the conventional matrix 
approach was extended to include internal finger reflections 
[3, 4]. The model divides an IDT into !4 wavelength 
sections. These sections are further divided into zones. 
Two of the zones are un-metalized areas (1/8 of a 
wavelength) around one zone that is comprised of a metal 
finger (1/4 of a wavelength). Each zone is modeled by a 
transmission line matrix equivalent circuit (Fig. 4). Two 
identical circuits model the un-metalized areas, while the 
middle circuit models the area under the metal finger. The 
transmission matrix relates the SAW the voltages Vj and 
V 2 , to the currents 1 1 to I 2 . The acoustic wave is assumed to 
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have entered from the left and travels through the element 
towards the right. In this model Z u and Z m are the acoustic 
impedances for the metalized and un-metalized areas, Co is 
the capacitance for a single finger, 9 U and 9 m are the 
acoustic angles of the substrate, and the turns ratio of the 
transformer is assumed to be 1 : 1 . 
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Where v m is the acoustic wave velocity under the 
metalized area and X is the wavelength of the device. The 
metalized velocity (v m ) is 3134 m/s for ST cut Quartz. 

The matrix (7) calculates the parameters for the 
metalized area, but cannot be used for the un-metalized 
sections. This leads to the transmission matrix (R u (f)) for 
the un-metalized region as is given by: 


[W)] = 


cosh(/Y/ (/)) Z u sinh(jO u (f)) 
^sinh(//?(/)) ash 


( 11 ) 


The transmission matrix (11) is determined by the 
acoustic angle 9 U and the un-metalized region’s acoustic 
impedance Z u . The acoustic impedance Z u is calculated 
with: 
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Fig. 4. Crossfield Model using Mason Equivalent Circuit for a !4 
wavelength section of an IDT. 


The transmission matrix for that represents the middle 
circuit of Fig. 4 for a metalized region that is assumed to be 
lossless is given by: 


K(/)]= 


cosh( AC/)) Z m sinhijdjf)) 
ys\rb(j0 m (f)) cosh(/<9 m (/)) 


(7) 


The transmission matrix (7) is determined by the 
acoustic angle 9 m and the metalized region’s acoustic 
impedance Z m . The acoustic impedance Z m is calculated 
with: 
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where k 2 is the piezoelectric coefficient, C s is the 
capacitance for a single pair of electrodes, H a is the aperture 
width, and f m is the frequency of the acoustic wave under 
the metalized area. The acoustic angle of the substrate 9 m , 
is given by 
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where f is the frequency, and f ra is the frequency of the 
acoustic wave under the metalized area: 


where fo is the synchronous frequency of the acoustic wave 
under for the un-metalized area. The acoustic angle of the 
substrate 9 U , is given by 
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where fo is the synchronous frequency of the acoustic wave, 
which is calculated using the acoustic wave velocity under 
the un-metalized area. The un-metalized velocity (v) is 
3158 m/s for ST cut Quartz. 

To find the transmission matrix for the Vi wavelength 
periodic element (RT(f)) one must multiply the three 
matrices together for both metalized region and the un- 
metalized regions adjacent to it: 

[*,(/)] = [*„(/}][*„(/}][*,(/)]• (14) 

To find the transmission matrix (Q(F)) for an entire IDT 
one simply raises the (Rt( 1)) matrix to the power of 2N p : 

[/Dr, (/)] = [*, (/)f\ (15) 

Where N p is the number of electrode pairs, so 2N p is the 
total number of electrodes in the IDT. 

The matrix for a SAW delay line is simply the 
multiplication of the matrices for the two IDTs and the 
delay or space between the IDTs. The SAW matrix is 
given by: 

[&4 W (/)] = [IDT, (/) D, (/) IDT, (/)] . 
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Fig. 5. Comparison of model results with data from a double side polish 
wafer, with 250 nm of aluminum. 

5 PROTOTYPE DESIGN 

A simple SAW delay line that consists of two identical 
un-apodized IDTs was chosen as a prototype to illustrate 
the capability of the models. Each IDT has 63 fingers that 
are 17 pm wide. The spacing between the fingers is 17 pm 
also. The center or synchronous frequency is 46.44 MHz, or 
a wavelength of 68 pm. The aperture height is 2730 pm. 
The delay length between the IDTs is 10 wavelengths or 
680 pm. The design was fabricated on two different quartz 
wafers. One with a single side polished, and one wafer 
with both sides polished. The aluminum thickness is 50 nm 
for the wafer with a single side polished and 250 nm for the 
wafer with both sides polished. 

6 RESULTS 

All three methods adequately model the frequency 
response amplitude for the main lobe and the first and 
second side lobes for cases without any mass loading due to 
the metal fingers (Fig. 6). The modified matrix more 
accurately captures the frequency shift due to the mass 
loading of the metal fingers. In Figure 5, the ideal first 
order model and the conventional matrix results are both 
centered about the synchronous frequency. The measured 
results and the Modified matrix results are both shifted 
down in frequency due to velocity changes from mass 
loading effects. 

When comparing the measured data from the two 
figures, it is noticeable that the main lobe peak of Fig. 6 
does not have the same artifacts as are seen on the peak of 
the main lobe in Fig. 5. These artifacts are due in part from 
bulk waves that are reflected from the polished bottom 
surface of the wafer. The roughness of the non-polished 
surface disperses the bulk waves which results in 
diminished artifacts in the main lobe peak (Fig. 6). Also 
note that the peak is not shifted as far in Fig. 6 as it is in 
Fig. 5. 



Fig. 6. Comparison of model results with data from a single side polish 
wafer, with 50 nm of aluminum. 

7 CONCLUSIONS 

In this paper, results from prototype devices are 
compared to the Impulse Response model, a conventional 
matrix model, and a modified matrix model. Ignoring mass 
loading from the metal fingers and the associated frequency 
shift, all three methods model the main lobe and next two 
side lobes fairly well. However, the results show that the 
modified matrix methods more accurately modeled second 
order effects such as frequency shift due to the metal 
thickness of the IDT fingers. Neither the Impulse response 
method nor the conventional matrix method model the 
frequency shift. Therefore, the modified matrix method is 
the most accurate for a wider range of parameters such as 
finger metal thickness. Future work will include a 
comparison the Coupling of Modes, the Impulse Response 
model, and the modified matrix model. 
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